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Abstract 

We construct a new algorithm for estimating the star formation rate (SFR) of galaxies from their 
infrared (IR) luminosity by developing the theory of the IR emission from a dusty H II region. The derived 
formula is SFR/(M® yr" 1 ) = {3.3 x 10- 10 (1 - r?)/(0.4 - 0.2/ + 0.6 e)}(L°^ s / L m ), where / is the fraction 
of ionizing photons absorbed by hydrogen, e is the efficiency of dust absorption for nonionizing photons 
from OB stars, and rj is the cirrus fraction of observed IR luminosity. The previous conversion formulae of 
SFR from the IR luminosity is applicable to only the case where the observed IR luminosity is nearly equal 
to the bolometric luminosity (starburst galaxies etc.), except for some empirical formulae. On the other 
hand, our theoretical SFR is applicable to galaxies even with a moderate star formation activity. That is, 
our simple and convenient formula is significantly useful for estimating the SFR of various morphologies 
and types of galaxies — from early elliptical to late spiral and irregular galaxies, or from active starburst 
to quiescent galaxies — as far as they have neither an abnormal dust-to-gas ratio nor an evident active 
galactic nucleus. 
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1. Introduction 

There are a number of studies estimating the star for- 
mation rate (SFR) of galaxies by using various observa- 
tional quantities. Indeed, we can estimate the SFR from 
colors of galaxies, hydrogen recombination lines, ultra- 
violet (UV) continuum and so on (Kennicutt 1998a). Of 
course, infrared (IR) radiation is also a very important 
tracer of the SFR. This is because some fraction of ioniz- 
ing and noninonizing photons from star forming regions 
are absorbed by dust grains which exist in or nearby 
these regions, and their energy is reradiated in IR spec- 
tral range (Soifer et al. 1987). 

In fact, most active star forming galaxies emit almost 
all of their radiative energy in IR (Soifer et al. 1987). For 
such galaxies, we can assume that their observed IR lumi- 
nosity is nearly equal to their bolometric luminosity. By 
adopting this assumption and a proper population syn- 
thesis model, Kennicutt (1998b) and Thronson & Telesco 
(1986) have estimated SFRs from the IR luminosity for 
starburst galaxies and for active dwarf galaxies, respec- 
tively. Likewise, taking into account the energy balance 
of such active star forming galaxies, we can also calculate 
their SFR from the IR luminosity (Devereux & Young 
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1991). 

For normal galaxies with a moderate star formation 
activity, however, the situation is more complex, because 
we can no longer assume their IR luminosity to be their 
bolometric one. In such cases, only empirical methods 
have been made to date. For example, Buat & Xu (1996) 
determined the SFR in terms of the far-IR luminosity, 
calculating the ratio of it to the UV luminosity. How- 
ever, their SFR is valid for late (Sb-) spiral galaxies only, 
since the ratio of far-IR to UV luminosity for early spi- 
ral (Sa-Sab) galaxies systematically differs from that for 
late spiral galaxies. Thus, we will establish more gen- 
eral method for the estimation of the SFR from the IR 
luminosity in this Letter. 

In general, we consider that massive young stars are 
formed in H II regions. In many H II regions, IR radia- 
tion is detected and it has a good correlation with ther- 
mal free-free radio emission, which comes directly from 
the ionized regions and is directly related to massive star 
formation (e.g., Spitzer 1978). Therefore, dust grains ex- 
ist in or nearby H II regions, and they absorb a part of 
radiation from young massive stars and reradiate this ab- 
sorbed energy in IR wavelengths. Petrosian et al. (1972) 
estimated the IR luminosity radiated from dusty H II re- 
gions with a simple analytic approximation. Our analysis 
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is based on their result. 

Unfortunately, not all the IR radiation originates from 
star forming regions in galaxies. The observed IR radi- 
ation from dust is the sum of the "warm" component 
distributed in or nearby star forming regions and the 
"cool" component (called cirrus component) distributed 
diffusely in the interstellar medium far from such regions 
(e.g., Helou 1986). Thus, the SFR calculated from the ob- 
served raw IR luminosity is clearly overestimated. When 
we estimate the SFR from IR luminosity, thus, the cir- 
rus component must be subtracted. In this Letter, we 
will subtract the cirrus component according to Lonsdale- 
Persson & Helou (1987), which presented a model of the 
cirrus fraction by using the ratio of IRAS 60/im and 
100/im fluxes. Also, we can estimate this fraction by 
applying a proper radiative transfer model for reproduc- 
ing multi- wavelength data of galaxies (e.g., Silva et al. 
1998, Efstathiou et al. 2000). However, the fraction of 
the cirrus component remains still uncertain. 

We stress again that only empirical studies on SFR 
from the IR luminosity of galaxies with a moderate star 
formation activity have been made to date. Thus, we de- 
velop the result of Petrosian et al. (1972) and construct a 
new algorithm to derive the SFR from the IR luminosity 
in this Letter. Thereby, we obtain a convenient formula, 
and calculate the SFR of galaxies from their IR lumi- 
nosity. This Letter contains the following sections: we 
review the IR luminosity of a dusty H II region in §2, our 
formula of SFR is constructed in §3, we discuss some re- 
sults derived from our formula in §4, and the conclusions 
of this Letter are summarized in the last section. 

2. Infrared Luminosity of H II regions 

First of all, let us look closely at the wavelength range 
of IR radiation from H II regions. We consider IR radi- 
ation to be mainly thermal black body radiation of dust 
heated by young massive stars in such regions. In fact, 
IR spectral energy distribution (SED) observed in H II 
regions can be fitted by the black body radiation of about 
30-50 K (Soifer et al. 1987). By the way, the temper- 
ature of dust within ionized regions is higher than that 
of dust out of there. That is, the peak wavelength of 
dust-IR radiation depends on the distance of dust grains 
from ionizing stars. It is, however, complex to take ac- 
count of more than two temperatures of dust and this is 
beyond the scope of this Letter. Instead of determining 
exact SED, we consider the dust luminosity in the whole 
range of 8-1000 /im, which covers almost all the emission 
from dust. The increment owing to stars at the shorter 
wavelength in this range is assumed to be subtracted. 

Next, we see the result of Petrosian et al. (1972). 
When the Case B approximation, which is an assumption 
of the large optical depth for every Lyman-emission-line 
photon (e.g., Osterbrock 1989), can be applied to a dusty 



H II region, Petrosian et al. (1972) derived the following 
equation for the dust-IR luminosity, Lf^ st (8 — 1000/im), 

Lf^(8 - 1000 M m) = L(Lya) + (1 - f)(hu) ion S 

+ £„o„ion(l-e- r ). (1) 

Here, L(Lya) is the luminosity of Lyman-a emission line, 
/ is the fraction of ionizing photons absorbed by hydro- 
gen, {hv)i on is the averaged energy of ionizing photons, 
S is the number of ionizing photons from central sources 
per unit time, £ n0 nion is the luminosity of nonionizing 
photons, and r is the optical depth of dust for nonioniz- 
ing photons. 

They assumed that all Lyman-a photons produced by 
hydrogen recombination processes in a H II region are 
absorbed by dust within this region and reemitted as IR 
radiation because of their large optical depth (i.e. Case 
B). This is reflected in the first term of the right hand side 
of equation (|]J). The second term represents the energy 
of ionizing photons absorbed directly by dust within the 
ionized region. Then, the last term denotes the energy of 
nonionizing radiation absorbed by dust in or nearby ion- 
ized region, especially, in molecular clouds surrounding 
this region. 

We now consider the luminosity of the Lyman-a emis- 
sion line within H II regions. Under Case B approxi- 
mation, every ionizing photon will eventually form one 
hydrogen atom with the n = 2 level. In this process, 
about two-thirds of recombining electrons will reach the 
2p state and go down to Is, emitting a Lyman-a pho- 
ton. The remaining one-third of recombining electrons 
will reach the 2s state, and two continuum photons will 
be emitted simultaneously within 1 second because the 
transition from the 2s to the Is is forbidden for any one 
photon process (e.g., Spitzer 1978). Thus, L(Lya) in 
terms of S is 0.67/i^Lya/<5', where ^Lya is the frequency 
at the Lyman-a emission. 

The luminosity of ionizing photons from central 
sources, L lon , is written as (hf}i on S. When we set 
hv^ya—10.2 eV and {hv)- lon ~ 15 eV, the luminosity of 
Lyman-a is given by 

L(Lya) = 0.67-J^L/L lon ~ 0.45/L ion . (2) 

Therefore, equation ([!]) is reduced to 

L^ st (8 - 1000/mi) = (1 - 0.55/)L ion + e£ nonion , (3) 

where e = 1 — e~ T is an averaged dust-absorption- 
cfficicncy of nonionizing photons from central sources in 
H II regions. 

3. Derivation of Star Formation Rate 

We will derive a new convenient formula to determine 
the SFR from the observed IR luminosity. Our result 
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is presented in equation ( |13[ ) in this section. First, in 
order to estimate the recent SFR, we assume that all 
of the total energy radiated from star forming regions 
is emitted by only young massive stars with the spectral 
types of O and B on the main-sequence. We fit the mass- 
luminosity relation of main-sequence stars from the table 
3.13 in Binney & Merrifield (1998) as 



log Z (771) 



log to + 4 for to > 20M m 
4 log m for to < 20M® 



(4) 



where m is the stellar mass in solar unit and /(to) is the 
stellar luminosity in solar unit as a function of m. 

We now must choose a initial mass function (IMF), 
tp(m), but the choice of it may affect our calculation of 
the SFR. There are two uncertainties in adopting a spe- 
cific IMF: the slope of ip(m), and the upper and lower 
cut-off of m. In this Letter, we adopt the Salpeter's 
IMF, ip(m) oc (to/M®) -2 - 35 (Salpeter 1955) and the mass 
range of 0.1 - 100 M®, since they seem to be widely ap- 
plicable. We normalize the IMF by J ip(m)dm = 1. We 
will discuss the uncertainty about the choice of an IMF 
more closely in §4. 

Next, we consider the mass range of O stars is 20 - 100 
iW® and that of B stars is 3 - 20 . Assuming that the 
total luminosity from star forming regions, L^ tai , is the 
bolometric luminosity of O and B stars, £qb, we obtain 

itotai = ios = « / l(m)^(m)dm, (5) 
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where a is a normalization so that M* = a J mip(m)dm, 
when we represent M* as the stellar mass formed newly 
in star forming regions. Since we adopt the Salpeter's 
IMF, the bolometric luminosity of O stars and that of B 



stars are determined as ^otype 



0.8L h o i and L™ pc 



0.2Lqb, respectively. 

Let us discuss what type of stars emits ionizing pho- 
tons dominantly. According to a model calculation by 
Panagia (1973), main-sequence stars of B type hardly 
emit these photons whereas those of O type radiate them. 
However, not all photons from O type stars have shorter 
wavelength than Lyman limit (912A). The energy frac- 
tion of ionizing photons from 05 stars is 0.6 and that of 
09 stars is 0.2 (Panagia 1973). Then, if we adopt 60 M 9 
and 20 Af® as the mass of 05 and 09 stars, respectively, 
and when the energy fraction of ionizing photons from 
stars of mass m is denoted by a(m), we get the following 
relation approximately, 



a(m/M®) = 0.01(to/M®). 



(6) 



Using equation (g), we can calculate the ionizing lumi- 
nosity of main-sequence O type stars. Then, we get 



T ion fi tf bol 

^Otype — U -°^Otypc 



0.4LQg. Now we assume that 



B type stars radiate nonionizing photons only. That is, 

T . _ r ion 
^lon — -^Otypei 



(L bo1 



Otype 



j ion \ _i_ r bol 
^OtypcJ T Li 



Btypc- 



Thus, Li on and £ non ion in equation (||) is reduced to 
0AL%& and 0.6L^, respectively. We now obtain the fol- 
lowing expression of the dust-IR luminosity of star form- 
ing regions, Lf£ st (8 - lOOO^m), 



1000/^m) = (0.4 - 0.2/ + 0.6 e)£gg. (7) 



In general, the observed far-IR luminosity contains 
both the component originated from star forming regions 
and that from the diffuse cirrus (e.g., Helou 1986). When 
we represent the cirrus fraction as to the total luminosity 
of star forming regions in a galaxy is, therefore, given by 



rSF _ rbol _ 
^total — -^OB — 



1-7? 



0.4 - 0.2/ + 0.6 e 



L°r(8- lOOOfim), (8) 



where L°* s (& - 1000/im) is the observed dust-IR lu- 
minosity of this galaxy. Furthermore, since M* = 
a J rmp(m)dm, the normalization a is removed by using 
equation (||): 



Mm Jg 1 ^ l(m)i>(m)dm As 



r bol j bol 

^ = 1.1x10-3^.(9) 



We now need to consider the time-scale of star for- 
mation because the SFR is the mass M* divided by this 
time-scale if the SFR is assumed to be constant during it. 
Since we are interested in the recent SFR of galaxies, we 
regard it as the representative life time of main-sequence 
OB stars. The stars on the main-sequence generate their 
radiative energy by the nuclear hydrogen reaction. Thus, 
the life-time of this phase is given by the following func- 
tion of mass of stars, 



t(m) = 0.0067 x 



f3mc 2 



(10) 



where c is the light speed, (3 denotes the fraction of 
mass consumed during this phase and Urn) is the mass- 
luminosity relation given by equation (^|). According to 
Schwarzchild (1958), the mass fraction, j3, through the 
CNO cycle is 0.13. Therefore, we obtain the life time of 
main-sequence as the following, 



t(m) m 
-—^ = 1.3 x 10 10 
yr 



M®J \l(m) 



(11) 



In this Letter, we choose the life time of star formation 
as the luminosity- weighted average life time of OB main- 
sequence stars. This is calculated by the following equa- 
tion, 



Jsm^ t{m)l{m)i>(m)dm 



(12) 



then, we obtain (t) = 3.3 x 10 6 yr. 
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Finally, from equations (|^) and (^), we derive the fol- 
lowing formula for the SFR of a galaxy in terms of its 
dust-IR luminosity by adopting 3.3 x 10 6 yr of the star 
formation time-scale, 

SFR _ 3.3 x 1CT 10 (1 - rj) Lj£ 8 (8 - 1000/an) 
M^yr- 1 _ 0.4 - 0.2/ + 0.6 e h$ '' 

This is the main result of this Letter. 

4. Discussions 

As stated in §1, the previous works on the SFR in 
IR luminosity is not applicable to galaxies with a mod- 
erate star formation activity, because we can no longer 
assume their IR luminosity to be their bolometric one 
for such galaxies. On the contrary, our equation ( |l3| ) 
is also reasonable for them. As a demonstration of our 
equation (|l3|), hence, we apply it to galaxies with a mod- 
erate star formation activity. We can adopt the suit- 
able sample with the data set of the IRAS luminosity, 
l fir( 40 - 120/mi), in Usui ct al. (1998), for example. 
Their sample contains 15 early (Sa-Sab) spiral galaxies 
showing rather high star forming activity in far-IR, and 
their averaged IRAS luminosity is about 5.8 x 10 9 L®. 

For absorbing fraction by neutral hydrogen, /, we 
adopt 0.26 derived by Petrosian et al. (1972) for Orion 
nebula. The efficiency of dust absorption for nonioniz- 
ing photons, e, is estimated to be 0.6 from the averaged 
1000-4000A extinction curve of the Galaxy (Savage & 
Mathis 1979) and the average visual extinction of Usui's 
sample (Ay = 1 mag from private communication with 
Usui). Here, we should note that when our formula is 
applied to the individual molecular clouds, these param- 
eters may depend on the geometry of clouds. Moreover, 
we choose 0.5 for the cirrus fraction, rj, of Usui's sam- 
ple, according to a model of Lonsdale-Persson & Hclou 
(1987). In addition, we have converted L^ s (8- 1000/zm) 
to 1,^(40- 120/zm) calculated easily from IRAS 60 and 
100 /im fluxes by the factor of 1.4, assuming the modi- 
fied black body radiation [ly oc vB Vl where l Vl i>, and 
B v are intensity, frequency, and the Plank function, re- 
spectively) of 30 K from an optically thin dust medium. 
Of course, this factor depends on the dust temperature. 
For example, if we assume 50 K and 15 K, the factors are 
changed to 1.1 and 5.6, respectively. Thus, we must note 
this point when our formula is applied to galaxies whose 
IR emission is dominated by almost only cold dust (e.g., 
15 K) owing to their quiescent star forming activity. 

Accordingly we obtain 2 yr _1 as the average value 
for the sample galaxies in Usui et al. (1998). This SFR 
is very similar to the value, 1.4 M®yr~ l , for the same 
sample via the equation for starburst galaxies in Ken- 
nicutt (1998b). Thus, we find coincidence between the 
two estimated SFRs. This coincidence originates from 



the following two reasons. [1] The coefficient of conver- 
sion from IR luminosity to SFR tends to decrease due 
to the cirrus fraction of IR. [2] It also tends to increase 
due to the lower optical depth of dust. Thus, we suggest 
that since the effect of the cirrus fraction is very effec- 
tive, that of the lower optical depth of dust is offset for 
adopted early spiral galaxies. 

Now we will discuss what type of galaxies our result 
can be applied to. Our algorithm starts from the ex- 
pected dust-IR luminosity of an H II region. There, we 
assume that the dust content is enough to absorb Lyman- 
a photons before they run away from this region. For 
most observed H II region, the Case B approximation is 
applicable, thus, all Lyman-a photons will be absorbed 
by dust within this region, even with a relatively low 
amount of dust (Spitzer 1978). Hence, equation ([l^) can 
be applied to galaxies with star formation in H II regions 
if they have neither an extremely small dust-to-gas ra- 
tio nor an active galactic nucleus (AGN). However, we 
can adapt equation (|l3|) to the AGN plus star forming 
galaxies, subtracting the AGN component of dust-IR lu- 
minosity from the observed IR luminosity. On the other 
hand, the limit of dust-to-gas ratio for the application of 
equation ( |l3| ) will be examined in our future work. 

Finally, we discuss the choice of an IMF, and an upper 
and lower cut-offs of mass. We examine the effect of vari- 
ance of the IMF, the upper limit mass or the lower limit 
mass. When we change only the IMF for the Scalo's IMF 
obtained from Binney & Merrifield (1998), the coefficient 
of the equation (|l3|) is smaller than that of the Salpeter's 
IMF by the factor of 0.8. If we substitute 60 M m into 
the upper limit mass, the recent SFR for the sample of 
Usui et al, (1998) becomes 1.4 times larger than that of 
100 M®. If the lower cut off mass is set 1 M®, the coef- 
ficient of the equation ([l3|) decrease by the factor of 0.4. 
In summary, the choices of an IMF, an upper limit mass 
and a lower limit mass cause uncertainties of the factor 
of about 2 to our result. 

5. Conclusions 

We summarize the conclusions reached in this Letter: 

1. Starting from the dust-IR luminosity expected the- 
oretically from a dusty H II region with the Case B ap- 
proximation, we formulate equation (|l3| ) to estimate the 
SFR of galaxies from their observed IR luminosity. 

2. This new simple formula for SFR contains three pa- 
rameters explicitly: [1] the fraction of ionizing photons 
from young massive stars absorbed by neutral hydrogen 
in star forming regions, /, [2] the averaged efficiency of 
dust absorption for nonionizing photons from young mas- 
sive stars, e, and [3] the cirrus fraction of observed IR 
luminosity, rj. 

3. Using equation ( |l3[) and adopting a proper set of 
three parameters, the recent SFR averaged for the sam- 
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pie in Usui et al. (1998), which consists of 15 early 
(Sa-Sab) spiral galaxies with moderate IR luminosity 
(~ lO 9 " 10 ^), is calculated to be about 2 M m yr" 1 . 

4. The derived convenient equation (|ll]) is applied to 
any galaxies forming young massive stars in H II regions 
as far as they have neither an extremely small dust-to- 
gas ratio nor an AGN, if we choose a set of parameters 
reasonable for applied galaxies. 

5. The coefficient of equation (|l^) has uncertainty of 
the factor of about 2 by the choices of a specific IMF, 
and its upper and lower limits of stellar mass. 
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